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Band structure is theoretically studied in partially flattened carbon nanotubes within an effective- 
mass scheme. Effects of inter-wall interactions are shown to be important in non-chiral nanotubes 
such as zigzag and armchair and can essentially be neglected in chiral nanotubes except in the close 
vicinity of non-chiral tubes. In fact, inter-wall interactions significantly modify states depending 
on relative displacement in the flattened region in non-chiral tubes and can convert semiconducting 
tubes into metallic and vice versa. They diminish rapidly when the chiral angle deviates from that 
of the zigzag or armchair tube, although the decay is slower in the vicinity of armchair tubes. 
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I. INTRODUCTION 

Carbon nanotubes were first found in a form of 
multi-wall cylinders, each of which consists of a rolled 
graphene sheetA single-wall nanotube, fabricated 
laterhas a unique electronic property that it changes 
critically from metallic to semiconducting depending on 
its tubular circumferential vector. This characteristic 
feature was first predicted by means of tight-binding 
models^ — and was successfully described in an effective- 
mass approximation^i^“— Experimentali^“— as well as 
computational studies^”— have discovered that large di¬ 
ameter nanotubes have an additional stable flattened 
structure. The purpose of this work is to study electronic 
structure of collapsed carbon nanotubes for arbitrary chi¬ 
rality within the effective-mass approximation. 

The observation of fully collapsed multi-wall car¬ 
bon nanotubes was reported in transmission electron 
microscopy^^— atomic force microscopyand scan¬ 
ning tunneling microscopy^^“— Multi-wall nanotubes 
was shown to exhibit structural deformations in FET 
devicesRecently, high-yield fabrication of high qual¬ 
ity collapsed tubes was reported, using solution-phase 
extraction of inner tubes from large-diameter multi-wall 
tubes^^ 

Actually, it is shown theoretically by first-principles 
energy minimization that both flattened and cylindrical 
nanotubes are stable or meta-stable and the energy of 
flattened tube is lower than cylindrical tubes with large 
diameterThe cylindrical nanotubes collapse into 
flattened tubes with a barbell-like cross section under 
hydrostatic pressure or in the presence of injected charge 
shown by molecular dynamics simulationsElec¬ 
tronic states were studied for collapsed armchair tubes in 
a tight-binding model^ and for collapsed zigzag tubes by 
density-functional calculationswhich demonstrated 
drastic modification in the energy region close to the 
Fermi level due to inter-wall interaction. 

Transport of crossed nanotube junctions results from 


interacting individual tubes and has been studied both 
experimentally^-— and theoretically^^-— The conduc¬ 
tance is found to depend strongly on the crossing an¬ 
gle with large maxima at commensurate stacking of lat¬ 
tices of two nanotnbes j^^i^^ A deformation of crossed 
carbon nanotubes, which may significantly affect the 
tunneling conductance between nanotubes, has been 
calculatedFurthermore, a pseudogap has been 
predicted to appear for an orientationally ordered crystal 
of nanotubes due to inter-tube transfer 

Effects of inter-wall interactions in multi-wall nan¬ 
otubes were also studied. In general, the lattice struc¬ 
ture of each nanotube is incommensurate with that of 
adjacent walls i^^i^^ This makes inter-wall electron hop¬ 
ping negligibly small as a result of the cancellation of 
inter-wall coupling in the absence of disorderIn fact, 
inter-wall hopping integrals vary quasi-periodically from 
site to site and their average over the distance of the 
order of the circumference vanishes. This property was 
extensively used for theoretical calculations of excitons in 
double-wall nanotubesFurther, it is closely related 
to very weak interlayer interactions in twisted bi- and/or 
multi-layer graphenes. 

Experimentally, each layer of some of epitaxially fab¬ 
ricated graphenes having many layers is known to be¬ 
have almost as a monolayerFurther, the elec¬ 
tronic structure of twisted bilayer graphene with nearly 
incommensurate lattice structure, both theoretically 
calculated^i^-— and experimentally observed 
shows a linear band dispersion near the charge neutral¬ 
ity point, suggesting weak interlayer interaction. On the 
contrary, the interlayer interaction drastically changes 
electronic states in displaced bilayer graphene having a 
commensurate lattice structure^iS^ The end of bilayer 
graphene can be closed and was observed experimen¬ 
tally after thermal treatment Geometry and elec¬ 
tronic structure of bilayer graphene with a closed edge 
were studied by a density functional calculation^iSi 

This paper is organized as follows: In Sect. [TTl an ef- 
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FIG. 1: (a) A schematic illustration of a collapsed carbon nanotubes and (b) its development map on graphene sheet. In (b), 
the coordinates system {x\y') and origin O’ are hxed onto the graphene sheet and the coordinates system {x^y) and origin O 
vary depending on the structure of a nanotube. The flattened region is denoted by the shaded area. 


fective potential of inter-wall interaction is derived in 
an effective-mass scheme. In Sect. im modification of 
band structure due to collapse is analyzed by perturba¬ 
tion of inter-wall interaction first for armchair and zigzag 
nanotube and its dependence on nanotube chirality and 
stacking in flattened region are discussed based on domi¬ 
nant terms. Numerical results are shown in Sect. lEI and 
a short summary is given in Sect. [Vl 


II. COLLAPSED CARBON NANOTUBES 

We consider a nanotube partially flattened as illus¬ 
trated in Fig. [T] (a). The width of the flattened region is 
denoted hy Lf/‘2^ and that of the curved region by Lcl‘2.. 
We have 

Lp + Lc = L, ( 1 ) 

where L is the circumference. Figure [ 1 ] (b) shows the 
development map. The tube is usually specified by chiral 
vector L, corresponding to the circumference, i.e., L = 
|L|. The direction of L measured from the horizontal 
direction is called the chiral angle and denoted by r]. 

In Fig. [T] (b), the right hand side of the line passing 
through the point O at (C cosr^, C sin?^) and perpendic¬ 
ular to L is folded down to form the lower half of the 
flattened nanotube. The coordinate of the point in the 
lower side of the flattened region corresponding to point 
r in the upper side will be denoted by f. Obviously, f 


is given by the mirror reflection with respect to the line 
perpendicular to L. In the nanotube, we shall use the 
coordinates (x, y) fixed onto the tube and therefore we 
have f = (—X, y) for r = (x, y). The coordinates of r and 
f in the coordinates (x', y') fixed onto the graphene sheet 
can be straightforwardly obtained. 

Figure [ 2 ] (a) shows the lattice structure of graphene, 
two primitive translation vectors a and b, and three vec¬ 
tors Ti (/ = 1,2,3) connecting nearest-neighbor atoms. A 
unit cell contains two carbon atoms denoted by A and B. 
In a tight-binding model, the wave function is written as 

tAW = ^ V’A(R)<^(r-R)+ ^ V’B(R)<A(r-R), (2) 

R,=R.a R,=FtB 

where 0(r) denotes a tt orbital. The amplitude -0 at 
atomic sites R = R^ or R^ satisfies 

3 

e'tpAi'R-A) = -7oT]V’b(Ra - n) 

+Y,V{RA,R'A)i’A{R'A) V(Ra,R'b)V’a(R'b), (3) 

R'a R'b 

3 

e-tpsiR-B) = 

y(RB, R'Ji^liR'^) + T/(Rb, R'b)V’b(R'b),(4) 

R'a R'b 

where 70 is the hopping integral between nearest- 
neighbor atoms within the wall and inter-wall hopping 
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FIG. 2: The lattice structure of graphene (a) and the first 
Brillouin zone (b). The primitive translation vectors are de¬ 
noted by a and b (|a| = |b| = a) and the vectors connecting 
nearest-neighbor atoms are denoted by fi, fb, and ts in (a). 
The reciprocal lattice vectors are denoted by a* and b*, and 
K = (27r/a)(l/3, 1/VS) and K' = (27r/a)(2/3, 0) in (b). 


integral T^(R, R') is nonzero only when carbon atoms at 
sites R and R' are very closely located in the opposite 
side of the flattened region. Since tt orbitals are symmet¬ 
ric within the wall, ^(R, R') is a function of |R — R'| 
well inside the flattened region. 

In a monolayer graphene the conduction and valence 
bands consisting of tt orbitals cross at K and K’ points of 
the Brillouin zone shown in Fig. [2] (b), where the Fermi 
level is locatedJ^iiS^ For states in the vicinity of the 
Fermi level 5 = 0, the amplitudes are written as 

(R^) + (5) 

^b(Rb) = (Rb) + '(Rb), 

( 6 ) 

with oj = exp(27ri/3)^ Envelope functions , F^, 
F^ , and Ff are assumed to be slowly varying in the 
scale of the lattice constant. The effective-mass approxi¬ 
mation is valid and well reproduces electronic properties 
as well as the band structure for energy range given by 
|e| < 37o>i^^— 

In the absence of inter-wall interactions, the envelope 
functions for the K point satisfy 

7i^(k)F^(r) =eF^(r), (7) 

= ‘o")' ® 

where k± = kx iky^ k = — iV, and 7 = \/3a7o/2 is 
the band parameter with lattice constant a = 0.246 nm. 
Here, wave vector k is measured from the K point de¬ 
noted by K. For the K’ point, we should exchange 
and k- in the Hamiltonian, i.e.. 


We shall construct a nanotube in such a way that the 
hexagon at L = n^a + n^b with integers and is 
rolled onto the origin. For translation r —> r + L, the 
Bloch function at the K and K’ points acquires the phase 


exp(zK • L) = exp ( -h ^ ), 

(11) 

exp(zK • L) = exp ( ^ ^, 

(12) 

where z/ = 0 or ±1, determined by 


ria rib = 3A/^ + u, 

(13) 


with integer N. Correspondingly, the boundary condi¬ 
tions for F(r) are given by 


F^(r + L)=exp(-^)F^(r), (14) 

F^'(r + L) =exp ^^^F^'(r), (15) 

Let T be the primitive lattice translation vector in the 
axis direction. 


T = maa + m6b, (16) 

with integers rria and m^. We have 

prria =na- 2n6, prrib = 2na - rib, (17) 

where p is the greatest common divisor of ria — 2nb and 
2na — rib. For translation r^r + T, the Bloch function 
at the K and K’ points acquires the phase 

exp(iK • T) = exp ( + ^ (ig) 

exp(iK' ■ T) = exp ( - ^ (19) 

where /x = 0 or ±1, determined by 

rria^mb = 3M -b /i, (20) 


with integer M. This shows that the K and K’ points 
are mapped onto k^ and k^ , respectively, with 



27r/i 


,K' ^ 

^ “ 3T ■ 


( 21 ) 


within the one-dimensional first Brillouin zone [—tt/T, 
+tt/T], with T = |T|. 

For the K point, the energies and corresponding wave 
functions are given hy^~— 


£naik) = F{k- (22) 

F^fes(r) = + i{k - , (23) 


where k is measured from the K’ point, i.e., K'. 
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with 




^ nks 


A = 


1 

71 V 

K^{n)-i{k-kA 

^K^inY + {k - kA"^ 


(24) 

(25) 

(26) 


where k is the wave vector in the axis direction, measured 
from the center of the one-dimensional Brillouin zone, n 
is an integer, s = +1 and —1 for the conduction and 
valence band, respectively, and A is the tube length. For 
the K’ point, we should replace u with —z^, /i with —/i, 
and k) with k) = ky. The band 

structure is illustrated in Fig. [3] in the vicinity of 5 = 0. 

In the presence of inter-wall coupling, the envelope 
functions satisfy the Schrodinger equation, 

?i(k)F(r) + J dr'V{r,r')F{r') = £F(r), (27) 


with 


Then, the effective potential of inter-wall interaction is 
explicitly given by 

VAAir)=J2 \[9{r-TiA)+9{r--RA]V{KA,-R-'A) 

Ra,R^ 

/ giK-(R^-RA) gi?7gi(K'-R^-K-RA) \ 

VabA) = E ^ [aA - ^a) + 9{r - R'B)]y(R^,R)j) 
Ra,R^ 

/ _^g*^g*K-(RB-R-A) gi(K'-R'^-K-RA) \ 

^ g-i»7e»K'-(R-B-R-A) )’ (^6) 

VbaA) = E ^ [5(r - Rb) + 9 {r - R'^)] V(Rb, R:4) 
Ra,R^ 

gi(K-R^-K'-RB) gir,giK'-(R^-RB) I’AO 

VBB(r)= E Ib(r-RB)+5(f-R'B)]V(RB,Ri5) 

Rb ,R^ 

/ giK-(R's-RB) _^-lg-ir7gi(K'-R'^-K-RB) \ 

^ ( _^gi77gi(K-R's-K'-RB) giK'-(R's-RB) ) * 



The (4,4) matrix effective-potential of inter-wall interac¬ 
tion is given by 

V(r,r') = ^(i’)'5(i’'- f), (30) 

with 



We should note that Ha and and also K and K' 
are in the coordinate system x'y' fixed onto the develop¬ 
ment map and r and f are in the coordinate system xy 
fixed onto carbon nanotubes. Thus, r and f should be 
converted into the x'y' system in the above equations. 

We have introduced a smoothing function gir) which 
varies smoothly in the range |r| < a and decays rapidly 
and vanishes for |r| ^ It should satisfy the condi¬ 
tions: 

E5(r-RA) = E^'(’'“^b) = 1, (39) 

Ra Rb 

J drg{r-RA) = J drgA-RB) = ^o, (40) 


where y^^(r), etc. are (2,2) matrices given by 






(32) 


etc. Explicit expressions for the effective potential are 
more easily written down in terms of 


E(r) 


fVAAir) EAB(r)\ 

v^ra(i*) vsBiy J ' 


(33) 


where is the area of a unit cell given by f2o = \/3a^/2. 
The function ^(r —R) can be replaced by a delta function 
when it is multiplied by a function such as F(r) varying 
smoothly in the scale of the lattice constant, i.e., ^(r — 
R) ^ — R)- 

The effective potential satisfies 


V{v) = V{vY, (41) 


and therefore 


where Vaa^)^ etc. are (2,2) matrices given by 


V:4A(r) 


tf/(r) VAa^'A)\ 
VIA A) VA;^'(r)J’ 


V(r, r')’*' = E(r)(5(r' — r) = V(r')S(r — r') = V(r', r), 

(42) 

which insures that V(r, r') is a Hermitian operator. The 
(34) Hamiltonian should satisfy the time reversal invariance 
under operation given b y^^'^d^^ 


etc. We can obtain V from E(r) by a simple unitary 
transformation. 
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FIG. 3: Schematic illustration of energy bands for (a) an armchair tube and (b) metallic and (c) semiconducting zigzag 

nanotubes with = 0. 


where (Tx is the Pauli spin matrix and ^ is an arbitrary 
phase factor. Thus, the inter-wall potential V (r) should 
satisfy 


0 CT;, 

CTz 0 


(Jx 0 
0 -CF. 


V{v) 

Viv)* 


0 cr^ 

Gz 0 


Gx 0 

0 -Gx 


= C(r), 

= C(r). 


(44) 

of the 


F(r) 

(45) 

/ ( 

trans- 

4 


lational properties: 


W^(r + L) = W^(r) 




V 


KK’ 


(r + L) =C^^'(r), 


V 


K'K 


(r + L) =r^'^(r), 


(46) 


V 


K'K' 


(r + L) = ^ (r) e 


K'K' / \ +47rizy/3 


and 


V^^{r + T) = V^^ir), 
V^^'iv + T) = v^^'iv) 
V^'^{v + T) = 

C^'^'(r + T) = V^'^\y). 


(47) 


Therefore, it can be expanded into a Fourier series, such 
that 


V 


KK 


2u\ 


(r) = ^ C„^^exp ^ (^4 - y - 


27rim 


-y 




L L 


27ri / 2/i\ 

—— [m - )y 

TV 3 


irK'Kf \ S^^rK'K \2Tlin 2711 
V (r) = 2 ^ Vn,m exp X 


L 

27ri , 


ZTTt ( ' 

- m H- ]y , 

TV 3 ) \ 

-r'rK' K' / \ X T^rK' K' [ZlTt ( 2v \ 27ri777/ 

^ (r) = ^ exp [— (^n + y jx + 


(48) 


with integers n and m. This shows that the inter-wall 
coupling gives rise to interactions among bands with 
same k value in the one-dimensional Brillouin zone. 

We shall expand the wave functions in terms of those 


1 


X, (49) 


ex.p[iK^ {n)x — ikjfy] 0 \ 

0 exj[)[iK,^ {n)x — ikjf y] j 

with Gm = 27rm/T. Then, we have 
7i[n, k^Gm]'Pnm^ E m'F n'm' = sFnm, (50) 


with 

H[n, k] = 
and 

^n-\-n' ,m—m' 


'H^[K^{n),k-k^] 




n+n ,m—m' n+n ,m—m' 


yK'K 

, n+n',m—m' 


t>yy 

n+n ,m—m' , 


(51) 


(52) 


For actual numerical calculations, the inter-wall hop¬ 
ping integral C(Ri,R 2 ) is chosen as^i2Sii22iii22i“Ji^ 


r(Ri,R2) = - 

|t|-6 


7i / 


[a^exp( 

<5 ) 


j(pi-t)(p2-t) 

7 oexp(^- ■^F^j[(pi-u)(p 2 -u) -H (pi-v)(p 2 -v)] , (53) 

where b is the distance between neighboring carbons in 
graphene, i.e., b = a/v^S, c the lattice constant along the 
c axis in graphite given by c/a = 2.72, and 6 the decay 
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FIG. 4: Some examples of the structure of the double-wall region of a flattened zigzag nanotube. The A and B sites remain 
the same, but the K and K’ points are exchanged between the upper (red) region and the lower (green) region, (a) — na 

corresponding to an AA stacked bilayer, (b) na < (^(n (c) C = (^ + |)<^? (d) C — (n ^)a corresponding to another AA 

stacked bilayer. 


rate of the tt orbital. Further, 71 is the hopping integral 
between nearest-neighbor sites of neighboring layers in 
graphite. Vectors, pi and p 2 are unit vectors directed 
along TT orbitals at Ri and at R 2 , respectively, t is a 
vector connecting the two sites, and u and v are unit 
vectors perpendicular to t and to each other. In the 
following numerical calculations, we use parameters 70 = 
2.7 eV, 71 = 0.4 eV, d/a = 0.185, and a = 1.4. 

The negative sign appearing in V (Ri, R 2 ) is due to the 
fact that the tt orbitals in the top and bottom sides of 
the nanotube have signs opposite to each other because 
of the tube geometry. Further, we choose the following 
smoothing function: 

j,(r) = ^e*p(-^). (54) 


III. WEAK INTER-WALL COUPLING 

A. Inter-Wall Potential in Zigzag and Armchair 

Figure m shows some examples of the structure of the 
flattened region in a collapsed zigzag tube (77 = 0). In 
zigzag tubes, A and B sublattices remain the same, but 
intra-valley components of inter-wall potential identically 
vanish, because the K point is mapped onto the K’ point 
and the K’ point onto the K point by the mirror reflection 
with respect to a line parallel to the axis. The effective 
potential in the flattened bilayer region becomes indepen¬ 
dent of position and is periodic as a function of C with 
period a/ 2 . 

We have in general 


with smoothing length d. This d is of the order of lattice 
constant a, but can be regarded as zero in the scale of 
the effective-mass approximation. In the following, we 
choose 2 < d/a < 5, for which the results are independent 
of this parameter. 

As shown by molecular dynamics simulations, coupling 
in the flattened region is slowly turned on in the vicin¬ 
ity of its edge 4^ When the coupling suddenly appears 
at a boundary, extra coupling terms may appear, being 
strongly localized at edges. In order to avoid such un¬ 
physical effects, we multiply the inter-wall hopping by 
the following function: 


0 (^X Xedge? ^edge) — ^ I 


^ X ^edge \ 


^edge ^ 

with erf(t) being the error function defined by 

2 


erf(t) = 


i: 


ds. 


(55) 


(56) 


The parameter Aedge describes the width of the region 
where coupling increases from zero to the value well in 
the flattened region. Actual calculations show that an 
extra effective potential localized at edges appears for 
Aedge^ 1 , but turns out to be negligibly small except 
in very narrow wires and therefore can safely be neglected 
for thick wires with collapsed structures. 


Vaa = —Vbb = -ui 



0 




-VbA = -U2 




(57) 


with real coefficients vi and V 2 and phases ipi and if 2 
varying with (/. Actually, difference ip 2 — is a rele¬ 
vant parameter changing the band structure, because a 
relative phase difference between the wave functions asso¬ 
ciated with the K and K’ point can be chosen arbitrarily 
and is not important. 

For the above potential, the band structure in the bi¬ 
layer region generally consists of two cone-like bands with 
crossing points displaced in the wave vector space and 
have different energies. These two cone-like bands repel 
each other when they cross. For C = ja/2 with j be¬ 
ing an integer, in particular, we have a bilayer with A A 
stacking and vi = 71 , V 2 = 0, and In this 

case, the two cone-like bands have different energies by 
± 7 i and do not interact each other. For C = ja/4, we 
have ip 2 = ^ 1 , for which two cone-like bands displaced 
from each other both in wave-vector and energy become 
independent. No AB stacked bilayer is formed. 

Figure [5] shows some examples of the structure of the 
flattened region for armchair tubes (77 = 7 r/ 6 ). Inter¬ 
valley components identically vanish, because the K and 













(a)2;=3n5 


(b) 3n/K2;<(3n+1/2)b (c) 2;=(3n+1/2)b 


(d)2;=(3n+i;b 


FIG. 5: Some examples of the structure of the double-wall region of a flattened armchair nanotube. The K and K’ points 
remain the same, but A and B sites are exchanged between the upper (red) region and the lower (green) region, (a) 2^ = Snb 
corresponding to AB stacked bilayer, (b) 3nb < 2( < (3n + |)^, (c) 2(^ = (3n + |)^ corresponding to AB stacking, (d) 
2(^ = (3n + l)b corresponding to AA stacking, with distance b = between neighboring carbon atoms within the plane. 


K’ points are mapped onto themselves after folding. An 
A site, however, turns into a B site and a B site turns 
into an A site, respectively. The effective potential is 
again independent of the position well inside the flattened 
region, but varies periodically as a function of ( with 
period 36/2, where 6 = a/\/3. We can set ( = 6(3j+p)/2, 
with integer j and 0<p<3. 

Numerical calculations show 


The matrix element 5'^/^ can be obtained from 
by replacing u with — z/. Further, and for 

different valleys can be obtained by setting = 0 in the 
above. 


B. Weak Inter-Wall Coupling: Zigzag Tube 


VaA = 


1 0 
0 1 


Vab = V^A = ^ 


VbB = —"^B 


1 0 
0 1 


1 0 
0 1 


(58) 


with real va^ vb^ and vab varying as a function of (. For 
p = 0 and 1, we have a bilayer with AB stacking and for 
p = 2, we have a bilayer with A A stacking. In fact, for 
p = 0 we have va = vab = 0 and = 7 i, for p = 1 we 
have Va = 7i and vb = vab = 0, and for p = 2 we have 
va = vb = 0 and vab = 

In the following, we shall consider effects of inter¬ 
wall interactions in the case of narrow flattened region 
Li?/L<Cl by perturbation analysis. This analysis is use¬ 
ful for understanding qualitative features of inter-wall in¬ 
teractions appearing in numerically obtained band struc¬ 
ture as shown in the next section. 

Because the effective inter-wall potential is indepen¬ 
dent of position, the spatial part of the matrix element is 
given by an overlapping integral. For KK elements, for 
example, we have 


nKK ^ ± 
^nn' 


(L+Lf)/4 ^(-L+Lf)/4 

, + / 

U (L-Lf)/4. J {-L-Lf)I4. j 

1 . fTT Lp 


^i[Ku{n)+Kr,{n')]x 

2 


Trn + n' — |z^ 


. [ttTf / ^ M 


X COS 


/ , 2 \ 
which for Lp/L^l becomes 

■( 


oKK ^ ^ 
^nn' 


L2 


n + n - -z/J 


(59) 


(60) 


For zigzag tubes, the effective inter-wall potential 
causes coupling between the K and K’ points. In the 
following, we shall consider the case that the flattened 
region has the structure of an AA-stacked bilayer, i.e., 
C = j {a j2) with integer j. 

In the case of semiconducting tubes {u = ±1), de¬ 
generate states associated with the K and K’ points are 
characterized by s' = 5 and \KF{n)\ = |/^_iy(n')|, giving 
n' = —n or i<i-y{n') = —Ku{n), as shown in Fig. [3] (c). 
The matrix elements are calculated as 

^vKK'] =[v^'Kn (61) 

L Jns,—ns L Jns,—ns T ^ 

independent of n. This shows that the two degenerate 
states split into two by the amount indepen¬ 

dent of bands. As will be shown in the next section, this 
can convert the tube into metallic for sufficiently large 

LfIL. 

In metallic case z^ = 0, there are two degenerate metal¬ 
lic linear bands for n = n' = 0 as shown in Fig. [3] (b). 
The matrix elements become 


KK'I 


J Os,Os 




i Os,Os' 


= 00 "(1 + ss') 




2L 


[y 

( 62 ) 

This shows that two degenerates states associated with 
the K and K’ points split into two by the amount 
± 7 i^f/( 2 I/) independent of bands and there is no band- 
gap opening. 

For parabolic bands n = ±no and n' = ±no with no > 0 
and 5 ' = s, the matrix elements are calculated up to 
linear order in /c, and the effective Hamiltonian within 
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the four degenerate states becomes 


-\-noK —rioK -\-noK' —noK' 

0 0 \ 

0 0 

0 0 

0 0 / 

(63) 

with S ~ {—l)^^k/Ko{no). With the use of the unitary 
matrix 




/ 

V 


1 

0 

0 


1 

0 

0 


0 0 \ 

0 0 / 


(64) 


the effective Hamiltonian is converted into 


= -71-^ 



0 

-i5 

0 

0 

-1 

0 

■i-iS 

+i6 

0 

+1 

0 

\ 0 

-i6 

0 

-1 


(65) 


This shows that the two bands, each doubly degenerate, 
are split by the amount ±71 Li? /L and then the remaining 
degeneracy is further lifted by S which is proportional to 
the wave vector and inversely to |n|. 


C. Weak Inter-Wall Coupling: Armchair Tube 


Armchair nanotubes have z/ = 0 and /i = ± 1 , and 
therefore are always metallic in the absence of inter-wall 
interaction, as shown in Fig. [3] (a). Thus, dominant inter¬ 
wall coupling is present only within each of the K and K’ 
points and . For p = 0 with AB-stacking 

structure, the matrix elements are calculated as 


L J ns,n s L J 


C 

-^7iS 


KK 
nn' • 


( 66 ) 


For n = n' = 0, there are two degenerate states s = ±1 
dit k = corresponding to the K point, and the effective 
Hamiltonian becomes 


'Weff = 






"fiLp 

2L 


2L 


2L 


-j\k - k^ 


"fiLp 

2L 


which gives the bands 


(67) 



This shows that the bottom of the conduction band re¬ 
mains at zero energy while the top of the valence band 
is lowered by jiLpjL. 

Because parabolic bands ±no with no > 0 and s' = s 
are degenerate, the effective Hamiltonian becomes 


= — 


")iLf 


(-l)no 1 \ 

1 (- 1)^0 ) ’ 


(69) 


giving energy shift of —[(—1)’^° ± 1]^iLf/{2L). This 
results in an alternate upward or downward shift by 
{—1)'^^^^^iLf/L for a band and no shift for another. 
For p = 1 with another AB-stacking structure, exactly 
the same results can be obtained by changing the phase 
of the wave functions in an appropriate manner. 

For the case of AA-stacking (p = 2 ) shown in Fig. [5] 
(d), the matrix element for the K point becomes 




KK^ 


s'[Ko(n) + i{k - kF\ 
s[Ko{n') - i{k - kf,)] 


Fi 


yjKo{n'y + {k- k^Y ' 


)5w • (70) 


For n' = —n including n = 0, we have Ko{n') = —Ko{n) 
and therefore. 


IvKK] , „/^ Ko{n) + i{k - k^) 'fiLp 

^ 2^ V«o(n)2 + (fc-fc^)2 L 

(71) 

For the metallic linear bands (n = 0), there are two de¬ 
generate states s = ± 1 . Off diagonal elements s' = —s in 
Eq. ([TT]) vanish, showing that there is no splitting. The 
bands are shifted by diagonal elements 


[V^^] 


Os,Os 


s{k - k^) jiLf 
\ k-k^\ L 


(72) 


This corresponds to a parallel shift in the negative k di¬ 
rection. In contrast, the linear bands at the K’ point shift 
in the positive k direction. 

For parabolic bands s' = s^ off-diagonal elements n' 
= —n 7 ^ 0 in Eq. (|7T]) causes the splitting of ±^iLf/L 
independent of the bands. Eor small k, diagonal elements 
for n' = n 7 ^ 0 become 




ns,ns 


\n 7iEf 

^ ^ |/^o(n)| L 


(73) 


corresponding to a parallel shift in k direction with 
amount decreasing with |n| as \n\~^. Sign of the par¬ 
allel shift is positive for odd n and negative for even n. 


D. Chiral Tubes: Dominant Terms 

It is known that essential properties of carbon nan¬ 
otubes can be specified by family index / given by 

f = 2na-nb. (74) 

Eor L having /, we have L • a = ^/a^, which means that 
the chiral vector lies on the line vertical to the horizontal 
axis with distance (l/2)/a from the origin. Erom Eq. 
Ci, we have na-\-ni) = Sua — /, meaning that the value 
of F is determined by /. Therefore, tubes with a given 
value of / have the same value of f and always metallic 
or semiconducting independent of individual values of Ua 
and rib. 
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FIG. 6: Some examples of the potential amplitude in the dominant-term approximation as the function of chiral angle rj with 
family index / = 144 defined in Eq. (EH), i.e., {ria^rib) = (72, 0), (73,2), ..., (96,48). The average of the absolute values of each 
element of V are shown for and ^ by (blue) solid lines and for and ^ by (red) dotted lines. The size of 

the one-dimensional unit cell T is also shown in units of L by (green) dashed lines. We chose relative displacement (/a = 0, 
1/4, and 1/a/3 = 0.5774- • • defined in Fig. [T] for three panels and a parameter Aedge is defined in Eq. dSSl). 


The dominant contribution of effects of inter-wall cou¬ 
pling in the flattened region may be estimated from the 
Fourier coefficients of small n and It is natural to 

choose n = 0 for both intra- and inter-valley terms. For 
intra-valley terms we choose m = 0 and for inter-valley 
terms we choose m = +/i for y^^ and m = —fi for 

yK'K^ 

In the following, in order to show the magnitude of 
the effective potential, we plot [L/^ etc. instead 
of ^ etc. themselves. Figure [ 6 ] shows the average of 

the absolute value of each element V separately for intra¬ 
valley (KK and K’K’) and inter-valley elements (KK’ and 
K’K) as a function of the chiral angle r] for tubes with 
family number / = 144, i.e., {na^nb) = (72,0), (73,2), 

• • - , (96,48). These tubes have = 0 and therefore are 
metallic. We have chosen the cases of C/a = 0, 1/4, and 
l/v^ = 0.5774 •••. 

The inter-valley coupling is significant only in the ex¬ 
treme vicinity of a zigzag tube 77 = 0. On the other 
hand, the intra-valley term gradually increases with 77 , 
with behavior strongly dependent on the width of the 
flattened region, L^jL. In fact, it oscillates with period 
roughly proportional to LjLp (T is a smooth and slowly 
increasing function of 77 ), but is not correlated with T 
that oscillates over wide range as shown in the figures. 
Further, we notice that the inter-wall interaction is es¬ 
sentially independent of relative displacement C except 
at 77 = 0 (zigzag) and 7 r /6 (armchair). This can be seen 
in the band structure itself as shown in the next section. 

Calculations are performed also for semiconducting 


tubes with family number / = 142, although not shown 
here. The behavior is qualitatively the same as in the 
case of / = 144, but the absolute value of the effective 
inter-wall coupling is smaller except in the case of 77 = 0 . 

In collapsed zigzag tubes with 77 = 0, inter-wall in¬ 
teraction is present only between the K and K’ points. 
This corresponds to the fact that phase factors, such 
as R^-k Ra) g 2 (K RA-K Ra) appearing in the 
off-diagonal elements of Vaa given in Eq. (j35j) and the 
corresponding terms in Vab^ etc. in Eqs. (|36|) - (|38|) . do 
not cancel out even after summation over R^, etc. 
When 77 slightly deviates from zero, however, these phase 
factors start to rapidly oscillate in a quasi-periodic man¬ 
ner because they involve K and K' considerably different 
(^ 27r/a) from each other. Thus, the effective inter-wall 
potential vanishes due to cancellation. 

In collapsed armchair tubes with 77 = 7 r/ 6 , on the other 
hand, the relevant phase factors involve same K or K'. 
When 77 slightly deviates from 7 r/ 6 , the phase factors start 
to oscillate in a quasi-periodic manner but the oscillation 
is relatively slowly-varying. Thus, the effective poten¬ 
tial remains nonzero for small Lp because of incomplete 
cancellation. Further, it decreases with the deviation of 
77 from 7r/6 more rapidly for wider flattened region, as 
shown in Fig. [6l 


IV. NUMERICAL RESULTS 

For actual calculations, we choose 71 as the energy 
unit. Further, we choose n = 0, ±1, • • •, ±7imax and for 
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FIG. 7: Calculated band structure of collapsed tubes with family index / = 144 having zigzag (a), its neighboring structures 
(b), armchair (f), and its neighboring structures (e) and (d). Parameters are defined in Eqs. ([T3l) . (ITgI) . and (|2D)), and illustrated 
in Fig. [U 
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FIG. 8: Calculated band structure of collapsed tubes with / = 144 having zigzag and its neighboring structures, (a) and (b) 
C/a = 1/4. (c) and (d) — 1/a/3- 


Tfl — 0, ihl, • • • , with TT-max — ^max — 10 ^ 15. 

This choice of the basis set gives convergent results at 
least for the bands lying in the zero-energy region in 
which we are interested. In the following, results for 
C/a = 0 , 1/4, and l/\/3 will be shown. Calculated band 
structure will be compared with that in the dominant- 
term approximation in which only the dominant term is 
taken as discussed in the previous section. 

Figure [3 shows examples of the band structure in the 
metallic case with / = 144 for C = 0. Figures [7] (a) 
and (b) present those in the vicinity of the zigzag struc¬ 
ture, i.e., {ua^ni)) = (72,0) and (73,2) corresponding to 
r]{7r/6)~^ = 0 and 0.046, respectively. The dominant- 
term approximation can describe the essential features 
of the bands near the Fermi level consisting of metal¬ 
lic linear bands, although wave vectors corresponding to 


zero energy are shifted and the velocity is lowered if we 
go beyond the dominant-term approximation. 

In zigzag nanotubes, the band structure is strongly 
modified by collapsing due to the strong inter-wall cou¬ 
plings, although the tube remains metallic because of the 
presence of linear bands at the Fermi level. In fact, the 
metallic linear bands associated with the K and K’ points 
are split in energy or shifted in the positive and negative 
k direction due to the inter-wall coupling as shown in 
Eq. (j62j) in the perturbation treatment in the previous 
section. The excited parabolic bands which are four-fold 
degenerate are first split into two sets at /c = 0 and then 
the remaining degeneracy is lifted by the k linear term 
as shown in Eq. (|65|). 

With the increase of 77 , i.e., when the structure de¬ 
viates from the zigzag case, effects of inter-wall in- 
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FIG. 9: Calculated band structure of collapsed zigzag tubes with / = 144 and 77 = 0. With the increase of Lf/F, the band 
structure gradually approaches that of an AA stacked bilayer with appropriately discretized wave-vectors in the circumference 
direction. 


teraction rapidly diminish. In fact, for = 

(73, 2), shown in Fig. [7] (b), the band structure is mod¬ 
ified due to the collapse in such a way that the effec¬ 
tive velocity in the axis direction is slightly reduced. 
This velocity reduction is in qualitative agreement with 
that observed experimentally^^i^Sri^iSZriffi and calcu¬ 
lated theoreticall y^^i^^~ — in twisted bilayer graphene. 
With the increase of 77 , however, the band rapidly be¬ 
comes unaffected by collapsing, although the results are 
not shown here. For chiral tubes, /i takes a nonzero value 
and therefore the K and K’ points become different in the 
one-dimensional Brillouin zone as denoted by short ver¬ 
tical dotted lines near zero energy. 

Figure [7] also shows results for tubes having a struc¬ 
ture close to T] = 7r/6 (armchair) and for = 0, i.e., 
(c) = (94,42), (d) (96,44), (e) (95,46), and (f) 


(96,48) corresponding to (c) 77 ( 7 r/ 6 )“^ = 0.893, (d) 0.930, 
(e) 0.965, and (f) 1. In the armchair tube, the figure 
shows results only for the K point and those for the K’ 
point are obtained by mirror reflection with respect to 
k = 0. 

In an armchair tube shown in Fig. [7] (f), the metallic 
band structure is strongly modified by collapsing and the 
tube becomes semiconducting due to band-gap opening. 
In fact, the bottom of the conduction band with n = 
0 remains at zero energy, while the top of the valence 
band with n = 0 is lowered roughly in proportion to 
LpjL due to inter-wall coupling. This is in qualitative 
agreement with Eq. (EH) obtained by the perturbation 
analysis. For excited parabolic bands, the qualitative 
features of effects of inter-wall coupling are in agreement 
with the perturbation analysis giving Eq. 
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FIG. 10: Calculated band structure of collapsed armchair tubes with / = 144 and r] = tt/G. With the increase of Lf/L, 
the band structure gradually approaches that of an AB stacked bilayer with appropriately discretized wave-vectors in the 
circumference direction. The tube becomes semiconducting due to inter-wall interaction in such a way that the bottom of the 
conduction band is fixed at zero energy, while the top of the valence band is lowered, forming a band gap. 


In chiral nanotubes effects of inter-wall interactions are 
considerably reduced and diminish with the decrease of 
T] from tt/G, although their decay is more gradual than 
in the vicinity of the zigzag tube. This has already 
been demonstrated in the behavior of the dominant terms 
shown in Fig.[6l The dominant-term approximation gives 
quite accurate results near zero energy, but starts to be¬ 
come less valid away from zero energy. 

The corresponding results for (/a = 1/4 and 1 /are 
shown in Fig. [8l Because the band structure is not af¬ 
fected by displacement ( in chiral nanotubes, only the 
results for zigzag and armchair nanotubes are shown. In 
zigzag and armchair nanotubes, the band structure de¬ 
pends significantly on (/a. 

As shown in Fig. SI the zigzag tube with r] = 0 has 


the structure of an AA stacked bilayer graphene in the 
flattened region for (/a = 0 and varies as a function of 
C with period a/2. For (/a = 1/4, the two layers are 
displaced from each other in a symmetric way, resulting 
in the symmetric band structure as shown in Fig. [8] (a). A 
small band gap appears for Lp/L = 1/4, but disappears 
for sufficiently large LpjL^ although explicit results are 
not shown here. 

The displacement C/a = l/v^S = 0.5774-•• in the 
zigzag case, shown in Fig. [8] (c), corresponds to the case 
that the top and bottom layers are slightly displaced from 
an A A stacked bilayer. This slight displacement results 
in repulsion between some bands of Fig. [7| (a), giving 
rise to the band-gap opening. This gap due to symmetry 
breaking is always present independent oiLpjL. 
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FIG. 11: Calculated band structure of collapsed tubes with / = 142 (semiconducting) having (a) zigzag and (b) its neighboring 
structure, and (c) and (d) near-armchair structure. The zigzag tube turns into metallic from semiconducting due to collapse 
as in (a). 


As shown in Fig. [5l for the armchair nanotube, the 
structure takes the form of AB stacking at C/a = 0 and 
C/a = l/(2\/3) = 0.2886 •••, and then the form of AA 
stacking Sit (/a = 1/\/3. This change repeats itself with 
period ajy/Z. Thus, in an armchair tube with Cja = 0.25, 
the structure is slightly displaced from the AB stacking. 
As shown in Fig. [8] (b), this slight displacement results 
in some distortion of the band structure of Fig. [7] (f) in 
such a way that the energy becomes asymmetric around 
the K point. In spite of the asymmetry, the tube remains 
semiconducting due to nonzero gap. 

For (/a = 1/a/ 3, the flattened region has the struc¬ 
ture of an A A stacked bilayer, and nanotubes become 
metallic independent of the width of the flattened re¬ 
gion, as shown in Fig. [8] (d), because linear bands cross 


the Fermi level. In agreement with Eq. (|7^. the metallic 
linear bands are shifted in the negative k direction. The 
parabolic bands are split and shifted in different k di¬ 
rections depending on band n qualitatively in agreement 
with Eq. (|73|) . 

In Eig.O the dependence on L^/A is shown for zigzag 
tubes with (/a = 0. With the increase of Lp/L^ the 
spectrum gradually takes a form of that of an AA stacked 
bilayer with appropriately discretized wave-vectors per¬ 
pendicular to the axis. The tube remains metallic in¬ 
dependent of LpjL. The dependence on the width of 
the flattened region for armchair nanotubes is shown in 
Eig. [To] for C/a = 0. The band structure again gradu¬ 
ally approaches that of an AB stacked bilayer. The band 
gap increases, takes a maximum, and then decreases with 
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LpjL, but always remains nonzero. 

As some examples for semiconducting nanotubes, we 
shall consider the case of / = 142 corresponding to 
{ua^nij) = (71,0), (72,2), ..., (94,46). In this case we 
always have u = —1 and /i = 0, i.e., the K and K’ points 
are both mapped onto the center of the one-dimensional 
Brillouin zone. Because qualitative feature of the de¬ 
pendence on the chiral angle is the same as in metallic 
nanotubes, we shall present results in the vicinity of the 
zigzag and armchair structure in Fig. [TTJ In fact, inter¬ 
wall interactions rapidly become small with the increase 
of T] from 7^ = 0. The same is true for i.e., inter¬ 

wall effects are most important for {ria^nh) = (94,46) for 
which T] = 0.977x(7r/6) closest to the armchair structure 
and decrease with the decrease of r] although more slowly. 

One most significant effect of the collapse is to convert 
semiconducting into metallic in the zigzag tube. This 
arises due to the splitting of two bands degenerate be¬ 
tween the K and K’ points due to inter-wall coupling as 
has been shown in the perturbation analysis, Eq. m- 
All chiral tubes remain semiconducting independent of 
Lf/L, although explicit results are not shown. 

V. SUMMARY AND CONCLUSION 

We have theoretically studied effects of inter-wall inter¬ 
action in collapsed carbon nanotubes within an effective- 
mass scheme. Inter-wall interactions in the flattened re¬ 
gion are represented by an effective potential connect¬ 
ing a point on the flattened region to its counter point. 
Effects of inter-wall interactions are most important in 
nonchiral nanotubes such as zigzag and armchair. In 
zigzag and armchair tubes, the band structure varies sen¬ 
sitively with the displacement, corresponding to the sen¬ 
sitive change of the band structure in bilayer graphene. 
In zigzag nanotubes, in particular, the collapsed tubes 
become metallic for sufficiently wide flattened region in¬ 
dependent of whether the uncollapsed tube is metallic or 
semiconducting. 


In chiral nanotubes, inter-wall interactions can essen¬ 
tially be neglected except in the close vicinity of zigzag 
and armchair tubes. Inter-wall interactions diminish 
rapidly when chiral angle deviates from 0 (zigzag) or 7r/6 
(armchair), although the decay is slower in the vicin¬ 
ity of the armchair tube. In fact, in chiral tubes closest 
to a zigzag and armchair tube, the semiconducting tube 
remains semiconducting even for very wide flattened re¬ 
gion. Such qualitative features of the chiral angle depen¬ 
dence can be understood through the magnitude of dom¬ 
inant terms corresponding to long-wavelength Eourier co¬ 
efficients of the effective inter-wall potential. 

Small band gap is inversely proportional to diameter 
in thick chiral semiconducting nanotubes, and is smaller 
than 7 i = 0.4 eV in nonchiral nanotubes, when band gap 
opens due to inter-layer interaction. Observation of these 
band gaps is required by means of precise measurement 
such as infrared transmission spectroscopy and scanning 
tunneling microscopy. Slightly reduced velocity due to 
the inter-layer interaction may be observed with scan¬ 
ning tunneling microscopy, angle-resolved photoemission 
spectroscopy, and Raman spectroscopy in the same way 
as in twisted bilayer graphene 

With the increase in the width of the flattened region, 
the band structure approaches that of a bilayer ribbon 
in which the electron motion in the ribbon-width di¬ 
rection is discretized under appropriate boundary con¬ 
ditions. Therefore, the band structure of collapsed nan¬ 
otubes can be obtained from a bilayer graphene by intro¬ 
ducing appropriate boundary conditions corresponding 
to the curved monolayer region. This problem is left for 
future study. 
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